Simultaneous ferromagnetic metal-semiconductor transition in electron-doped EuO 
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We present a general framework to describe the simultaneous para-to-ferromagnetic and 
semiconductor-to- metal transition in electron-doped EuO. The theory correctly describes detailed 
experimental features of the conductivity and of the magnetization, in particular the doping depen- 
dence of the Curie temperature. The existence of correlation- induced local moments on the impurity 
sites is essential for this description. 
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At room temperature stoichiometric europiumoxide 
(EuO) is a paramagnetic semiconductor which undergoes 
a ferromagnetic (FM) transition at the Curie tempera- 
ture of Tq — 69 K. Upon electron doping, either by O 
defects or by Gd impurities, this phase transition turns 
into a simultaneous ferromagnetic and semiconductor- 
metal (SM) transition with nearly 100 % of the itiner- 
ant charge carriers polarized and a sharp resistivity drop 
of 8 to 13 orders of magnitude, depending on sample 
quality [l|, 0, 0, 0|- Concomitant with this transition 
is a huge colossal magnetoresistance (CMR) effect @, 
much larger than in the intensely studied manganates 
Q . These extreme properties make electron-doped EuO 
interesting for spintronics applications. Known since the 
1970s, these features have therefore recently stimulated 
more systematic experimental studies with modern tech- 
niques and improved sample quality 0, 0> HI as well as 
theoretical calculations d, [l(| • 

In pure EuO the FM ordering is driven by the Heisen- 
berg exchange coupling between the localized Eu 4/ mo- 
ments with spin Sf = 7/2 [11]. Upon electron doping, 
above Tc, the extra electrons are bound in defect lev- 
els situated in the semiconducting gap, and the transi- 
tion to a FM metal occurs when the majority states of 
the spin-split conduction band shift downward to overlap 
with the defect levels. Although this scenario is widely 
accepted, several questions of fundamental as well as ap- 
plicational relevance have remained poorly understood. 
(1) Why does the magnetic ordering of the Eu 4/ sys- 
tem occur simultaneously [4( with the SM transition of 
the conduction electron system? (2) What is the order 
of the transition? While the magnetic ordering of the 4/ 
system should clearly be of 2nd order, the metallic tran- 
sition requires a finite shift of the conduction band and, 
hence, seems to favor a 1st order transition. (3) How 
can the critical temperature Tc be enhanced by doping 
for spintronics applications? While in the Eu-rich com- 
pound EuO i-z a systematic Tc increase due to the O 
defects (i.e. missing O atoms) is not observed experi- 
mentally P, [H, a minute Gd doping concentration sig- 
nificantly enhances Tc @, II3]- An O defect in EuOi-z 



essentially binds the two excess electrons from the extra 
Eu 6s orbital and, therefore, should not carry a magnetic 
moment. As shown theoretically in Ref. [lfj, the pres- 
ence of O defects with two-fold electron occupancy does 
not enhance Tc, in agreement with experiments [l], 0]. 
In the present work we focus on the Gd-doped system 
Eui_. y Gdy and calculate the temperature and doping de- 
pendent magnetization and resistivity from a microscopic 
model. We find that the key feature for obtaining a Tc 
enhancement is that the impurities not only donate elec- 
trons but also carry a local magnetic moment in the para- 
magnetic phase. 

The model. — A Gd atom substituted for Eu does not 
alter the Sf = 7/2 local moment in the Eu Heisenberg 
lattice but donates one dopant electron, which in the in- 
sulating high-temperature phase is bound in the Gd 5d 
level located in the gap. Therefore, the Gd impurities 
are Anderson impurities with a local level Ed below the 
chemical potential \i and a strong on-site Coulomb re- 
pulsion U > n — Ed which restricts their electron occu- 
pation essentially to one. The hybridization V with the 
conduction band is taken to be site-diagonal because of 
the localized Gd 5d orbitals. The Hamiltonian for the 
Eui-yGdyO system then reads, 
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where the first term in Eq. Jl} denotes conduction elec- 
trons with spin a. The Eu 4/ moments Si on the 
lattice sites i — 1,...,N are described in terms of a 
Heisenberg model H c f with FM nearest and next-nearest 
neighbor couplings Jij and an exchange coupling J c f to 
the conduction electron spin operators at site i, <?, = 
(V 2 ) Haa' c L^' c ia'> witn = Sk exp(ikx ; ) c kcr and 
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f aa i the vector of Pauli matrices. The Gd impurities at 
the random positions i = 1, Ni are described by H cd . 
For the numerical evaluations we take U — > oo for sim- 
plicity. 

For the present purpose of understanding the general 
form of the magnetization m(T) and the systematic dop- 
ing dependence of Tc it is sufficient to treat the 4/ 
Heisenberg lattice, H c f, on mean field level, although 
recent studies have shown that Coulomb correlations in 
the conduction band can soften the spin wave spectrum in 
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similar systems j!3lll4|. The effect of the latter on m(T) 



can be absorbed in the effective mean field coupling of 
the 4/ system, J4 / = ^ . . We therefore choose J4 / 
such that forpure EuO it yields the experimental value of 
Tc = 69 K [l|, 0, 0, 0|. For simplicity, we don't consider 
a direct coupling Jdf between the 4/ and the impurity 
spins, since this would essentially renormalize J c f only. 
The indirect RKKY coupling will also be neglected, since 
for the small conduction band fillings relevant here it is 
FM, like J^ , but much smaller than J^-. 

In the evaluations we use a semi-elliptical bare conduc- 
tion band density of states (DOS) with a half width Do = 
8eV, (consistent with experiment Q), centered around 
A ~ 1.05 D above the (bare) defect level Ed- The other 
parameters are taken as J4/ = ^ . = 7 • 10~ 5 Z?o> 
J cf = 0.05A), E d = -0.4L> , and V = ttV 2 = 0.05D 2 , 
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where J c f 3> J4/ because J4 / involves a non-local matrix 
element. 

Self consistent theory. — The averaging over the ran- 
dom defect positions is done within the single-site T- 
matrix approximation, sufficient for dilute impurities. 
This yields for the retarded conduction electron Green's 
function G ccr (k, uj) in terms of its selfenergy £ 0(T (u>), 

G ca (k,uj) = [w + /i-e k -E C(T (w)] _1 (4) 
Z ca (u J )=n I \V\ 2 G dr7 (u J )- J cf (S)a (5) 

where C?<j CT (u;) is the defect electron propagator and (S) 
the average 4/-moment per site. In mean field theory it 
is obtained, together with the conduction electron mag- 
netization m, as 
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where f{oj) is the Fermi distribution function and 
A ca -(uj) = — J^ k ImG cff (fc, u>)/tt the conduction electron 
DOS of the interacting system. In order to treat the 
strongly correlated spin and charge dynamics of the An- 
derson impurities without double occupancy beyond the 
static approximation, we use a slave particle representa- 
tion and employ the non-crossing approximation (NCA) 
[lit . For EuO the DOS at the Fermi level is so low or 
even vanishing that the Kondo temperature is well below 
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FIG. 1: (Color online) Conduction (upper panel) and impu- 
rity (lower panel) electron DOS per lattice site for T <C Tc- 
The impurity concentration is ni — 1.2% . The metallic phase 
is fully spin polarized. Inset: c-electron DOS in the param- 
agnetic phase for m = 1.2%, T = 100.8 K and for m = 0.6%, 
T — 88.0 K. The chemical potential lies in the gap. 

Tc and Kondo physics plays no role. In this high-energy 
regime the NCA has been shown to give quantitatively 
reliable results 16]. This remains true even for a finite 
magnetization, where the NCA would develop spurious 
potential scattering singularities near Tk only [17] . One 
obtains the following set of equations for Gd CT (w) in terms 
of the auxiliary fermion and boson propagators Gj a , Gb, 
their spectral functions Af a , Ab and their selfenergies 



S/a(w) = r / de [1 - /(e)] A ca ( £ )G b (uj - e) 
E h H = T Y,J dsf(e)A ca (e)G fa (u + e) 
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Gda{u) = J ^[G f<T (Lo + s)Ab(e)-A fa (e)Gl(e-Lu)] 

(10) 

Note that in Eqs. ([U A ca (e) is the interacting DOS, 
renormalized by the dilute concentration of Anderson im- 
purities and the 4/-spins according to Eq. (J3]) . For de- 
tails of the NCA and its evaluation see [16( . The equa- 
tions (Til fTUl) form a closed set of selfconsistent integral 
equations. They are solved iteratively, fixing the total 
electron number per lattice site in the system, 

CO (w) + ni Ada{w)\ = ni (11) 

by the chemical potential fi in each step. 

Electrical conductivity. — The current operator j can 
be derived from the continuity equation, dpi/dt + V -j = 
0, and the Heisenberg equation of motion for the to- 
tal local charge operator pi at site i. Because the 
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FIG. 2: (Color online) Simultaneous FM and semiconductor- 
metal transition as seen in the T-dependence of the magne- 
tization m and of the resistivity p — l/a. The Curie tem- 
perature depends significantly on the impurity concentration 
ni . The inset shows a typical experimental p(T) curve, taken 
from Ref. 0|. 

impurity Hamiltonians H c f, H^f conserve pi, only c— 
electrons contribute to the current, and one obtains [HI, 
j = (e/H) J2ka de^/Ok cj CCT c k(T . The linear response con- 
ductivity then reads for a local selfenergy (ljj . 
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Results and discussion. - The results of the self- 
consistent theory, Eqs. (|4T [Tl"j) , and for the conductivity, 
Eq. (|12[) . are presented in Figs.[TH3] They allow to draw a 
complete picture of the FM semiconductor-metal transi- 
tion in Gd-doped EuO. The spectral densities per lattice 
site above and below the transition are shown in Fig. [T] 
In the paramagnetic, insulating phase the hybridization 
between d- and c-electrons necessarily implies the ap- 
pearance of a conduction electron sideband (Fig. [TJ in- 
set), situated below fj, and at the same energies inside 
the semiconducting gap as the impurity d-band. The ri- 
band (not shown) has a similar width and shape as the 
c-sideband. The combined weight of the c-sideband and 
the d-band adjusts itself selfconsistently such that it just 
accommodates the total electron number, n = n.j. Note 
that the weight of the <i-band per impurity and spin is 
< 1/2, because the doubly occupied weight is shifted to 
C/^oo0. 

The c-4/ exchange coupling J c f induces an effective 
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FIG. 3: Conductivity a for various temperatures (left panel) 
and Curie temperature Tc (right panel) as function of impu- 
rity concentration m. The data points at m — in the left 
panel are extrapolations. 

FM coupling between the electrons of the c-d system. 
Hence, either the 4/- or the c-d-electron system can 
drive a FM transition, depending on which of the (cou- 
pled) subsystems has the higher Tc ■ We have chosen J c f 
(see above) large enough that the transition is driven by 
the c-d-electrons, because this will yield detailed agree- 
ment with the experiments 0, HIHj]- I n this case, Tc is 
naturally expected to increase with the impurity density 
nj. The results for the T-dependent conduction elec- 
tron magnetization m(T), Eq. ([7|), and for the doping 
dependence of Tc are shown in Fig. [2j lower panel, and 
in Fig. [3J right panel, respectively. It is seen that not 
only Tc increases with the impurity concentration, in 
agreement with recent measurements on Eui-yGdyOi-a, 
[3 EH' but also that m(T) has a dome-like tail near T c , 
before it increases to large values deep inside the FM 
phase. From our theory this feature is traced back to 
the mean-field-like 2nd order FM transition of the elec- 
tron system, while the large dome in the magnetization 
further below Tc is induced by the FM ordering of the 
4/ system, whose magnetization is controlled by J4 / and 
sets in at lower T. This distinct feature is again in agree- 
ment with the experimental findings and lends sig- 
nificant support for the present model for Eui_. y Gd u O. 



We note that the Eu-rich EuOi_ K samples of Ref. [12 1 
also show a magnetization tail and a Tc enhancement, 
suggesting (small) magnetic moments on the O defects. 
However, the nature of the O defects requires further ex- 
perimental and theoretical studies. The conduction elec- 
tron polarization P(T) = m(T) / n c (T) does not show this 
double-dome structure and below Tc increases steeply to 
P = 1 (not shown in Fig. The FM phase is connected 
with a spin splitting of the c- as well as the (i-densitics 
of states, as shown in Fig. [1] The narrow d-band in- 
duces a Fano dip structure in the c majority band and 
a small sideband in the c minority band. Note that for 
the present scenario the existence of preformed local mo- 
ments on the impurities, induced by strong Coulomb re- 
pulsion U , is essential. Without these moments the tran- 
sition of the electron system would be purely Stoner-like, 
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FIG. 4: (Color online) Conductivity ( — ) and magnetization 

( ) as a function of ju for T = 89.6 K and T = 99.2 K. 

At the ungated points, fi « -0.414 D for T = 89.6 K and 
fi « —0.410 Do for T = 99.2 K, the total electron number per 
lattice site is equal to the impurity concentration, n = rij. 

and, because of the extremely low conduction electron 
DOS at the Fermi level, its Tc would be far below the 
Curie temperature of the 4/ system, so that no doping 
dependence would be expected [10J. 

We now discuss the conductivity and the simultaneity 
of the FM and the SM transitions. In the paramagnetic 
phase, the system is weakly semiconducting, because \x 
lies in the gap (Fig.Q] inset). When the FM transition oc- 
curs, the impurity d-band must acquire a spin splitting in 
such a way that at least part of the minority ci-spectral 
weight lies above the chemical potential fi, in order to 
provide a finite magnetization. Since near the transition 
the spin splitting is small, the majority <i-band must, 
therefore, also be shifted to have overlap with \i (Fig. [TJ, 
and so must the hybridization-induced c-electron side- 
band (which eventually merges with the main conduction 
band for T sufficiently below Tc). This immediately im- 
plies a transition to a metallic state, simultaneous with 
the FM transition, as seen in Fig.[2j Because of the small, 
but finite thermal occupation of the states around (j,, we 
find that this shifting of spectral weight occurs continu- 
ously, which implies the FM semiconductor-metal tran- 
sition to be of 2nd order (see Fig. [3]). The doping m 
dependence of the conductivity is shown in Fig, [3l left 
panel. It is seen that the metallic transition can be driven 
by increasing if T > Tc{ni = 0). 

As an alternative to Gd-doping the charge carrier con- 
centration n can be controlled independently of the im- 
purity concentration nj by varying the chemical potential 
fx, e.g. by applying a gate voltage to an EuO thin film. 
The conductivity a and magnetization to as a function of 
\x are shown in Fig.[4]for two temperatures. To both sides 
of the ungated system (n = nj) a increases exponentially 
upon changing /i, characteristic for semiconducting be- 
havior. By increasing /i, the FM-metallic transition is 
finally reached. I.e. the magnetization can be switched, 
in principle, by a gate voltage. The non-monotonic be- 
havior of a towards more negative [i reflects the energy 
dependence of the c sideband. A more detailed study will 
be presented elsewhere. 



To conclude, our theory indicates that in Gd-doped 
EuO the existence of preformed local moments on the 
impurity levels inside the semicondicting gap is essential 
for understanding the distinct shape of the magnetiza- 
tion m(T) near the ferromagnetic semiconductor- metal 
transition. The FM ordering is driven by these impu- 
rity moments which are superexchange coupled via the 
4/ moments of the underlying Eu lattice. This scenario 
immediately implies an increase of the Curie tempera- 
ture with the impurity concentration, in agreement with 
experiments. The double-dome shape of m(T) arises be- 
cause of the successive ordering of the dilute impurity 
and of the dense Eu 4/ systems, as T is lowered. The 
dynamical accumulation of conduction spectral weight at 
the chemical potential, induced by the hybridization V 
and the constraint of an emerging magnetization at the 
FM transition, implies the FM and the SM transition to 
be simultaneous and of 2nd order. The magnetization 
can be switched by applying a gate voltage. This might 
be relevant for spintronics applications. 
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